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ABSTRACT 

Weibel instability turns out to be the a ubiquitous phenomenon in High-Energy Den- 
sity environments, ranging from astrophysical sources, e.g., gamma-ray bursts, to labo- 
ratory experiments involving laser-produced plasmas. Relativistic particles (electrons) 
radiate in the Weibel-produced magnetic fields in the Jitter regime. Conventionally, 
in this regime, the particle deflections are considered to be smaller than the relativis- 
tic beaming angle of I/7 (7 being the Lorentz factor of an emitting particle) and the 
particle distribution is assumed to be isotropic. This is a relatively idealized situation 
as far as lab experiments are concerned. We relax the assumption of the isotropy of 
radiating particle distribution and present the extension of the jitter theory amenable 
for comparisons with experimental data. 

Subject headings: radiation production; Weibel instability; laboratory astrophysics; 
high-energy-density physics; gamma-ray bursts; shock waves 



1. Introduction 



In general, high Mach number shocks, e.g., relativistic shocks of gamma-ray bursts , must be 
highly turbulent. It has been shown tha t the Weibel instability IWeibell (119591 ): I Fried! (|1959l ) is 
responsible for the GRB shock formation (jMedvedev &: Loeblll999i ) on the microscopic level. This 
instability is driven by the anisotropy of the particle distribution function (PDF) associated with a 
large number of particles reflected from the shock into the upstream region. This theore tical predic 



tion has recently been confirmed in a number of state-of-the-art numeric al sumulations (jSilva. et al 



20031 : iFrederiksen. et al.1 12004 : iNishikawa. et al.1 120031 : ISpitkovskvH2008r ). The Weibel instability is 
also observed in laser-plasma experiments. In particu lar, it is the goal of the H ercules experiment 



at th e university of Michigan ([Huntington et al.l 120081 : for technical details, see iMaksimchuk et al 



2008) to create and diagnose the Weibel instability and turbulence in the laboratory high-energy 
density plasmas, as a part of the Laboratory Astrophysics and High-Energy Density Physics pro- 
grams. 

The state of the Weibel turbulence corresponds to the self-organized nonlinear regime of the 
Weibel instability, which is characterized by reorganization of the currents and magnetic fields 
via strong interaction and hierarchical merger process of current filaments. The Weibel-generated 
magnetic fields are very small-scale, of order several plasma skin depths, c/u p , which is much 
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smaller than the typical Larmor radius of particles in such fields. Thus, if relativistic electrons are 
present in such a plasma, they will produce emission that is different from the standard synchrotron 
radiation. Such radiation, referred to as the "jitter radiation" has distinct spectral properties 
(|Medvedevl 120001 ) . has been predicted to be emitted from the Weibel turbulence. It has also been 
predicted that jitter radiation can explain various observational data from GRBs. At last but not 
least, jitter radiation has been suggested as an interesting diagnostic of the Weibel turbulence in 
laser plasma experiments. Here we elaborate more on jitter radiation theory in application to the 
experiments. 



Weibel turbulence 



The instability under consideration was first predicted by IWeibell (|1959l ) for a non-relativistic 
plasma with a n anis otropic distribution function, and the physical interpretation was provided 
later by iFriedl (|1959l ). who considered the extreme case of anisotropy — two counter-streaming 
particle (plasma) beams. In essence, the two electron-proton plasma streams experience deflections 
in seed magnetic fields due to the Lorentz force, e(v x B)/c, so that protons (and electrons) moving 
in opposite directions concentrate in spatially separated current filaments. The magnetic field of 
these filaments appears to increase the initial magnetic field fluctuation. The growth rate and the 
wavenumber of the fastest growing mode (which, in fact, sets the spatial correlation scale of the 
produced field) are of order the plasma frequency, u p , and the plasma skin depth, c/oj p , respectively. 
The current in the filaments and the associated magnetic fields increase until the energy density in 
the fields reaches about e B ~ 10% of the kinetic energy density of the streaming particles, which is 
enough to rapidly isotropize the particle distribution; hence the instability quenches. 

At longer times, the plasma with the Weibel-generated currents and fields enters the turbulent 
state, referred to as the "Weibel turbulence", in which current filaments begin to interact with each 
other, forcing like currents to ap proach each other and merge. The filament coalescence is a hier- 
archical and self-similar process (jMedvedev. et al.ll2005l ). For filaments with the initial separation 
~ Dq the magnetic field correlation length in the non-relativistic and relativistic filaments regimes 
are 

\ B (t) = D 2 t l { - 2T0 - NR \ X B (t) ~ ct. (1) 

Here, the typical non-relativistic time-scale is determined by Eq. (TjQ). The coalescence time may 
be written as 

( C / V ) 1n 4,-l m 

T0,NR ~= —p= ~ 10 W , (2) 

Here v is the be am velocity and we ass umed the typical value: e B ~ 10 _1 . Numerical PIC sim- 
ulations showed (jMedvedev. et al.ll2005l ) that both a non-power-law non-relativistic regime and a 
power-law regime are clearly present in the dynamics. The power-law fits yield X B (t) oc t a with 
a ~ 0.8. It should be noted that the field scale growth is somewhat analogous to the inverse cascade 
in two-dimensional magnetohydrodynamic (MHD) turbulence, with the crucial difference that the 
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former is an entirely kinetic process since at such small scales ~ c/lo p the MHD approximation is 
completely inapplicable. 



3. Jitter radiation 

The jitter regime realizes when the deflection angle, a, of a particle in the magnetic field is 
smaller than the relativistic radiation beaming angle ~ I/7, that is a -C A9. In this case, the 
velocity v of a particle is almost constant whereas its acceleration w = v varies with time. Calcu- 
lating the Fourier component of the electric field using the Lienard-Wiechart (retarded) potentials, 
one arrives at the following expression for the total energy emitted per unit solid angle df2 per unit 
frequency du>: 



dw e /UJ 



2irc 3 



n x 



n I x w (J 

c, 



dn—, (3) 



where = J we luJ>t dt is the Fourier component of the particle's acceleration, u' = u (1 — n • v/c), 
and n is the unit vector pointing towards the observer. We need to express the temporal Fourier 
component of the Lorentz acceleration, w = (e/'ymc)v x B, taken along the pa rticle trajectory i n 



terms of the Fourier component of the field in the spatial and temporal domains [MedvedeyJ ((2006). 
In the static case, i.e., when the magnetic field is independent of time, the ensemble-averages 
acceleration spectrum reads: 

(|w w f) = (2TTV)- 1 y |wk| 2 <5(c/ + k-v)dk, (4) 
|w k | 2 = (ev /~fmc) 2 (5 a/ 3 — v~ 2 v a vp) VK a ^(k), (5) 

K a p(r,t) = T _1 y _1 J B a (r' ,t')Bp(r' + r, t' + t) dr'dt is the second-order correlation tensor of the 
magnetic field. 

We adopt the following geometry: the Weibel current filaments are aligned with the z direction 
and their magnetic fields lie predominantly in the x—y plane. For a GRB, this geometry corresponds 
to a shock which is located in the x-y-plane and is propagating along z-direction. As the shock is 
propagating through a medium, the produced field is transported downstream (in the shock frame) 
whereas new field is continuously generated at the shock front. Thus, the field is also random in 
the parallel (z) direction. Similar structure of the magnetic fields is expected in the laboratory 
experiments with the field being random in the plane perpendicular to the beam propagation 
direction. Thus, Weibel turbulence shall be highly anisotropic. Both the theoretical considerations 
and realistic 3D simulations of relativistic shocks indicate that the dynamics of the Weibel magnetic 
fields in the shock plane and along the normal to it are decoupled. Hence, the Fourier spectra of 
the field in the x — y plane and in z direction are independent. Thus, for the Weibel fields at shocks, 
the correlation tensor has the form 



Kapfc) = C(5 afS - s a sp)f z (k\\)f xy (k_L), 



(6) 
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where s is the unit vector along the filaments (and normal to the shock front in the GRB case), C is 
the normalization constant proportional to (B 2 ), f z and f xy are the magnetic field spectra along s 
and in the perpendicular plane, respectively, k± = (k 2 + k 2 ) 1 / 2 and fey = k z , and finally, the tensor 
(Sap — $aSp) is symmetric and its product with s is zero, implying orthogonality of s and B. 



Numerical simulations (IFrederiksen. et al.l 12004 ) also indicate that the field transverse spec- 
trum, f X y, is well described by a broken power-law with the break scale comparable to the skin 
depth. We expect that the spectrum f z , has similar properties. Therefore, we use the following 
models: 

k 2ai fc 2a2 
/z(fc||)= {n\ + k 2 )^ fxy{k±)= (Kl+kl)to' (?) 

where K\\ and n± are parameters determining the location of the peaks in the spectra, ax, «2> P2 
are power-law exponents below and above a spectral peak {fix > 0,2 + 1/2 and > 012 + 1, for 
convergence at high- A;). 

We now evaluate Eqs. (J3J) , (JSJ) - The scalar product of the two tensors is 

(<W - v a vp/v 2 ){8 a p - s a sp) = 1 + {s a v a ) 2 /v 2 = 1 + cos 2 6, (8) 

where we used that 5 aa = 3. Here is the angle between the filament direction (and normal to 
the shock for GRBs) and the particle velocity (in an observer's frame), which is approximately the 
direction toward an observer, that is v||k for an ultra-relativistic particle (because of relativistic 
beaming, the emitted radiation is localized within a narrow cone of angle ~ 1/7)- Eq. (j4]) becomes 

(l w -'| 2 ) = ^(l + cos 2 6) jf z {k\\)f xy (k L )5{J + \L-v)dk\\d 2 k L . (9) 

Equations ([3]), ([9]) fully determine the spectrum of jitter radiation from relativistic electrons prop- 
agating through the Weibel turbulence. 



3.1. Isotropic electron distribution: the GRB case 

First, we can simplify the vector expression in ([3|). Indeed, in the ultrarelativistic case, the lon- 
gitudinal component of the acceleration is small compared to the transverse component, w\\/w± ~ 
I/7 2 -C 1. Therefore v and w are approximately perpendicular to each other. Second, the dominant 
contribution to the integral over d£l comes from small angles 9 ~ I/7 with respect to the parti- 
cle's velocity. Therefore, we approximately write u' ~ u (l — v/c + 6 2 /2j ~ ^uj (l — v 2 /c 2 + 9 2 ) = 
TjU) (6 2 + 7 -2 ). We now can replace integration over the solid angle dVt ~ 9 d9 dcf) with integra- 
tion over d(j)duj'/u) and integrate equation ([3|) over the azimuthal angle, eft, from to 2tt. The 
angle-averaged spectral power emitted by a relativistic particle moving through small-scale ran- 
dom magnetic fields, under the assumption that the deflection angle is negligible and the particle 



- 5 - 



trajectory is a straight line, has been derived elsewhere (jLandau k, Lifshitall97ll ; lMedvedevll2000l ): 



dW 



e 2 uj 
2^3 



W,„/ 



a/7 2 



+ 



uJ- 



duj'. 



(10) 



A spectrum from a shock viewed at an arbitrary angle, < < tt/2, is illustrated in Figure 
[H which represents full numerical solutions of Eqs. (|10p . ([7]), Q for three different viewing angles. 
In calculation of dW/ckv, the emitting electrons were assumed monoenergetic, for simplicity. An 
important fact to note is that the jitter radiation spectrum varies with the viewing angle. When 
filaments are (and a shock velocity is, in the GRB case) along the line of sight, the low-energy 
spectrum is hard F u oc v 1 , harder than the "synchrotron line of death" {F v oc u 1 ' 3 ). As the viewing 
angle increases, the spectrum softens, and when the filaments are orthogonal to the line of sight, 
it becomes F v oc v°. Another interesting feature is that at oblique angles, the spectrum does not 
soften simultaneously at all frequencies. Instead, there appears a smooth spectral break, which 
position depends on 0. The spectrum approaches ~ v below the break and is harder above it. 



3.2. Beam electron distribution: the lab case 

Unlike a GRB shock, the distribution of radiating electrons is anisotropic in m ost of lab experi 



merits. In particular, to diagnose the Weibel turbulence in the Hercules experiment ([Huntington et al 



2008), it has been suggested to launch a probe, nearly monoenergetic electron beam through plasma 
with the Weibel fields. Thus, the geometrical shape and the electron energy-momentum distribution 
of the probe beam are important. Here we assume the electrons to be monoenergetic and neglect 
the geometrical divergence of the beam, for simplicity (in the experiment, it is likely somewhat 
smaller than the relativistic beaming cone of I/7 anyway). 

We again start from the Lienard-Wiechart potentials and the expression for the emitted power 
([3]). We adopt the geometry such that the unit vector s is along the filaments, v is the particle's 
velocity and the unit vector n is toward an observed. We also define the unit vector v = v/v 
and (3. Unlike the isotropic case, we do not neglect terms proportional to |n • w^/j 2 although they 
are small compared to those proportional to (w^l 2 . Using a similar approach, we arrive at the 
following expression for the emitted power per frequency, per solid angle, per electron: 

dWu e 2 1 



dwdU (2tt) 2 c 3 (1 - p( n ■ v)) 4 

x P 2 ^J f z (k\\)f xy (k±)6 (u (1 - /3(n • v )) - k • v) dk\\d 2 h± 

x {(n-v) 2 (l + (s-v) 2 )+ [(s-n) 2 + (s-v) 2 -2(s-n)(s-v)(n-v)]}. (11) 



One can readily see that the emitted power along the electron probe beam, i.e., when n • v = 1, 



- 6 - 



is proportional to the electron acceleration spectrum: 

dW L 



duo dft 



oc Iw^r, (12) 



which is represented by the integral over k, up to a constant. Thus, the electron beam can directly 
probe and diagnose the structure of the magnetic field distributions: f z if the beam is aligned with 
the Weibel filaments and f xy if the filaments are probed "edge-on"; see Figure[T]for the w-spectrum. 
The angular shape of the radiation pattern (i.e., the angular distribution of intensity) can also be 
readily calculated by integrating over the frequency: 

^oc(l-/?(n.v))- 5 «(l + ( 7 ^) 2 )" 5 . (13) 

The latter expression is valid for small angles i? between the probe beam direction and the line of 
sight. 



4. Conclusions 

The jitter spectra can deliver much information on the structure of the Weibel magnetic fields. 
As one can see, the spectrum depends on the spatial spectra of the magnetic fields modeled by 
Eqs. (J7J). In particular, when viewing angles are and 90 degrees, the contributions of the parallel 
and transverse magnetic field spectra are decoupled. For instance, for 9 = 0, the peak of the jitter 
radiation spectrum and its high-energy asymptotic slope are uniquely determined by the parallel 
correlation length kh and the large-/c magnetic field spectrum slope with r\ = 2a\ — 2(3\. Similarly, 
the transverse jitter spectrum (at 9 = vr/2) allows one to deduce these parameters for the transverse 
magnetic field spectrum, f xy . At intermediate angles, one can determine the relative orientation of 
the current (and magnetic) filaments in the target and the radiation detector. 

It seems feasible to obtain jitter radiation in a laser-plasma experiment, such as Hercules. 
The Weibel turbulence to be studied will have much in common with the upstream region of a 
gamma-ray burst collisionless relativistic shock. It may be so even up to and at the main shock 
compression, where the Weibel filaments are present. In the downstream regio n the filaments ar e 



destroyed and the fields are significantly isotropized, as follows from simulations (|Spitkovskvl l2008) . 
In such a turbulence state, jitter radiation can still be present, but it will produce more isotropic, 
relatively soft spectra resembling those at 9 ~ ir/2. 

This work has been supported by NASA grants NNX07AJ50G, NNX08AL39G, NSF grant 
AST-0708213, and DOE grant DE-FG02-07ER54940. 
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Fig. 1. — The log — log plots of (w^ 2 vs k (thin lines) and of dW/dio vs cu (thick lines), for three 
viewing angles = 0, n/10, n/2. The axes units are arbitrary. In this calculation we used f z = f xy 
with a = 2, /3 = 20, k = 10, u = 1. The exponent C = C( a ^P) is model dependent. We also chose 
7 = 1 in order to align the peaks of (w^l 2 and dW/dio. Note that the actual peaks are at values 
k, <jj lower than 10 by a factor two or three. Note also that the spectrum dW/doj levels off at 
oblique angles at frequencies much smaller than Ki^sinG, whereas |wk| 2 indeed starts to flatten 
at k ~ kv sin O. 



